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Abstract 

When modelling metapopulation dynamics, the influence of a sin- 
gle patch on the metapopulation depends on the number of individuals 
in the patch. Since there is usually no obvious natural upper limit on 
the number of individuals in a patch, this leads to systems in which 
there are countably infinitely many possible types of entity. Analo- 
gous considerations apply in the transmission of parasitic diseases. In 
this paper, we prove central limit theorems for quite general systems 
of this kind, together with bounds on the rate of convergence in an 
appropriately chosen weighted l\ norm. 
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1 Introduction 



Metapopulations, introduced by Levins (1969), are used to describe the evo- 
lution of the population of a species in a fragmented habitat. The metapop- 
ulation consists of a number of distinct patches, together with (a summary 
of) the population present in each patch, and its development over time is 
governed by specified within and between patch dynamics. In the Markovian 
structured mean- field metapopulation model of Arrigoni (2003), the state of 
the system consists simply of the numbers of individuals in each patch. In- 
dividuals reproduce within patches and migrate between patches, and each 
patch is subject to random catastrophes, which reduce its population to zero. 
Letting N be the total number of patches, thought of as being large, and let- 
ting X^' 1 denote the number of patches with i individuals at time t, the 
transitions out of state X e Z+ + to states X + J in her model are as follows: 

J = e^-V - eW at rate Nix^d* + 7(1 - p)), % > 2; 

J = e (°)-e (1) atrateiVx 1 (di + 7(l-p) + /t); 

J = e( i+1 ) - e w at rate Nix%, % > 1 

J = - e w at rate Nx*k, % > 2 

J = e( fe+1 ) - e( fc ) + - e« at rate Nix i x k p-f, k > 0, % > 1 

here, x := N~ Y X. The total number iV := ^2j >0 X N, i of patches remains 
constant throughout, and the number of individuals in any one patch changes 
by at most one at each transition. The per capita death and birth rates 
(di), {pi) within each patch are allowed to vary with the current population 
size i, but in the same way in all patches; they would usually be chosen 
to correspond to one of the traditional single species demographic models. 
The per capita migration rate 7 is also the same for all individuals, as is the 
probability p that a migration is successful, and a successful migrant chooses 
its new patch uniformly at random. Each patch is independently subject to 
catastrophes at the same rate k. 

If there were an absolute upper limit for the number of individuals in each 
patch, the model would be a finite dimensional Markov population process. 
The behaviour of these finite dimensional models can be approximated using 
the methods pioneered by Kurtz (1970, 1971), who was able to establish a 
law of large numbers approximation, in the form of a system of ordinary dif- 
ferential equations, and a corresponding diffusion approximation. However, 
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there are no upper limits on population number in the usual single popula- 
tion models, and it is the stochastic evolution according to the rules of the 
model that dictates the region in which population numbers typically lie. 
Thus it seems unnatural to introduce an a priori upper limit in the system 
above, just because more than one population is being modelled. The same 
considerations surface in a number of other population models, including the 
epidemic models of Luchsinger (2001a,b) and Kretzschmar (1993), and the 
model of cell behaviour as a function of the copy number of a particular gene 
in Kimmel & Axelrod (2002, Chapter 7). Instead, it makes sense to consider 
Markov population processes with a countably infinite number of dimensions 
as models in their own right. 

A law of large numbers in a general setting of this kind was first estab- 
lished by Eibeck & Wagner (2003). Under appropriate conditions, Barbour 
& Luczak (2008, 2011) strengthened the law of large numbers by provid- 
ing an error bound, in a weighted t\ norm, that is close to optimal order 
in N. In this paper, these latter results are complemented by a central limit 
approximation, together with a corresponding error estimate. 

Our general setting, as in Barbour & Luczak (2011) [BL], is that of fam- 
ilies of Markov population processes X N := (X^, t > 0), N > 1, taking 
values in the countable space X + := {X G Z+ + ; J2 m >o^ m < °°}- The com - 
ponent X^^ of Xj* represents the number of individuals of type j that are 
present at time t, and there are countably many types possible; however, at 
any given time, there are only finitely many individuals in the system. The 
process evolves as a Markov process with state-dependent transitions 

X ->■ X + J at rate NajiN^X), X G X+, JeJ, (1.1) 

where each jump is of bounded influence, in the sense that 

J C jx G Z z+ ; ^2 \X m \ < J* < oo j, for some fixed J* < oo, (1.2) 

m>0 

so that the number of individuals affected at each transition is uniformly 
bounded. Density dependence is reflected in the fact that the arguments 
of the functions aj are counts normalised by the 'typical size' N. Writing 
1Z := R+ + , the functions ay. 1Z — > R + are assumed to satisfy 

< £ e ^o, (1-3) 

Jej 
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where 1Z := {£ 6 1Z: £j = for all but finitely many i}; this assumption 
implies that the processes X N are indeed pure jump processes, at least for 
some non-zero length of time. To prevent the paths leaving X + , we also 
assume that J 1 > — 1 for each I, and that = if = for any J E J 

such that J 1 = —1. 

In the finite dimensional case, the law of large numbers is expressed in 
terms of the system of deterministic equations 



| = E'Mf). (1.4) 



In [BL], it is assumed that 



£jaj(0 = A£ + F(£), (1.5) 

where A is a constant Z + x Z + matrix, and ( 11.41) is then treated as a perturbed 
linear system (Pazy 1983, Chapter 6). Under suitable assumptions on A, 
there exists a measure \i on Z+, defining a weighted £i norm || • on 1Z, and 
a strongly || • | ^-continuous semigroup {R(t), t > 0} of transition matrices 
having pointwise derivative R'(0) = A. If F is locally || • || M -Lipschitz, the 
solution x of the integral equation 

x t = R{t)x + [ R(t- s)F(x s )ds, (1.6) 
Jo 

for \\xo\\^ < oo, replaces that of ( 11. 4ft as an approximation to := iV -1 ^^. 
Under suitable conditions, it is shown in [BL, Theorem 4.7] that 



sup = 0(N-V 2 y/te£N), 

0<t<T 

except on an event of probability of order 0(iV _1 log N), provided that 
\\xq — a^olU = 0(N~ l l 2 ^\ogN). The conditions under which this approxi- 
mation holds can be divided into three categories: growth conditions on the 
transition rates, so that the a priori bounds, which have the character of 
moment bounds, can be established; conditions on the matrix A, sufficient 
to limit the growth of the semigroup R, and (together with the properties 
of F) to determine the weights defining the metric in which the approxima- 
tion is to be carried out; and conditions on the initial state of the system. 
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The conditions are described in the next section. They are all needed in the 
current paper, too, in which we investigate the difference x? — x t in greater 
detail. 

Our main result, Theorem 16. 1[ shows that, under some extra conditions, 
it is possible to construct a diffusion process Y on the same probability space 
as in such a way that 

sup ||iV 1/2 (xf -x t )-Y t \l = 0(iV- fcl ), (1.7) 

0<t<T 

except on an event of probability of order 0(N~ bi ), for specific values of 
h\ and 62- With the best possible control of moments, as for the model 
of Arrigoni (2003) mentioned above, one can take any b\ < 1/4 and any 
b 2 < 1, provided that the initial conditions are appropriately chosen. The 
process Y can be interpreted as the infinite dimensional analogue of the 
diffusion approximation in Kurtz (1971), satisfying the formal stochastic dif- 
ferential equation 

dY t = {A + DF(x t )}Y t dt + dW t . (1.8) 

Here, dW is a time-inhomogeneous white noise process with infinitesimal co- 
variance matrix cr 2 (t) := Xljej" JJ T aj(xt), and Y has time-inhomogeneous 
linear drift with coefficient matrix A + DF(x t ). In particular, if x is an equi- 
librium of the deterministic equations, satisfying Ax + F{x) = 0, then Y is an 
infinite dimensional Ornstein-Uhlenbeck process, with constant drift coeffi- 
cient matrix A+DF(x) and infinitesimal covariance matrix Xljejr JJ T &j{x)- 

Basic approach 

The structure of the argument is as follows. It is shown in [BL, (4.8)] that, 
under suitable conditions, the process x N satisfies an equation very similar 
to CLD: 

x? = R(t)x£+ [ R(t-s)F(x^)ds + m?, (1.9) 
Jo 

where 

t rt 

R(t-s)dmf = mf + / R(t - s) Am? ds, (1.10) 

Jo 
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and 

mf := x?-x»- /W + F(xf)}ds (1.11) 
Jo 

is a local martingale. Taking the difference between (jl.9p and (II. 6p . and 
multiplying by y/~N, gives 

C/f = R{t)U? + f R{t-s)DF{x B )[U?]ds + rj? + iV^mf, (1.12) 

with L/f := N 1/2 {xf - x t } and 

r/f := A^ 1 / 2 /* R(t — s){F(x^) — F(x s ) — DF(x s )[N~ 1 ^ 2 U^}} ds. (1.13) 
Jo 

Starting with this representation of U N , the first step is to show that rj N is 
uniformly small with high probability, so that the randomness in U N is driven 
principally by the process N l / 2 m N . This quantity is in turn determined, 
through (ll.lOp . by the local martingale N l ^ 2 m N . The next step is to show 
that N l l 2 m N is close to a diffusion W, formally expressible as 

W t := ^JWj(Aj(t)), (1.14) 
Jej 

where the {Wj, J G J} are independent standard Brownian motions, and 
Aj(t) := f Q aj(x s )ds: this is the diffusion appearing in (ll.8p . Analogously 

to (ll.lOp . we then show that we can define a process W such that 

W t := W t + f R(t-s)AW s ds, (1.15) 
Jo 

and that W is close to N l l 2 m N . Finally, returning to (I1.12p . we show that 
U N is close to the solution Y to the analogous equation 

Y t = R{t)Y + [ R(t-s)DF(x s )[Y s ]ds + W t , (1.16) 
Jo 

which in turn can be shown to exist and be unique. The random process solv- 
ing ( I1.16P at first sight seems rather mysterious. However, partial integration 
represents W t as J Q * R(t — s) dW s , and so the expression for Y can indeed be 



6 



interpreted as the variation of constants representation of the solution to the 
formal stochastic differential equation (II. 8p . 

In the remaining sections, this programme is carried out in detail. Sec- 
tion [2] is concerned with specifying the conditions under which the main 
theorem is true, and with recalling the results from [BL] that are needed 
here. In the subsequent sections, the steps sketched above are examined in 
turn. 

2 Assumptions and preliminaries 

We assume henceforth that (11. 2p and (jl.3p are satisfied. Since the index j G 
Z + is symbolic in nature, we fix an v e TZ, such that reflects in some 
sense the 'size' of j: 

> 1 for all j > and lim = oo. (2-1) 

We then assume that most indices are large and that most transitions involve 
some large indices, in the sense that, for 7m := {j: u(J) < M} and Jm '■= 
{J: Ji = for all j Tm}, we have 

\T M \ ^n^; \J M \<n 2 M^, (2.2) 

for some n\,ri2 and 0i < $2] note that in fact j3 2 < 1fi\J* also. As a 
consequence of these assumptions, for any s > j5\, there exists K s < oo such 
that 

X>0T < oo; J2^' S < K S M~ S+ ^- (2.3) 
moreover, if v(J) := max^j. jj^o} v {j)i then, for any s > (3 2 , 

J2»( j r s < E < j y s < km-**, (2-4) 

for some K' s < oo. 

Moment assumptions 

In the proofs that follow, it is important to be able to show that x N is 
largely concentrated on indices j with not too large. This is shown to 
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be the case in [BL, Section 2], under the following 'moment' assumptions. 
Defining S r (x) := ^> x j {u(j)} r , x G TZq, and then 



u r (x) ■.= X)«jKE^w)} r ); v ^ ■= E^( x )(E Ji Wi)y 



x ETZ, the assumptions that we need are as follows. 



(2.5) 



Assumption 2.1 For v as above, assume that there exist r^lx, ?max > 1 



such that, for all X G X + , 
JeJ 



i>o 



< oo, 



< r < r 



(i) 



(2.6) 



and also that, for some non-negative constants k r i, the inequalities 



U (x) < k 01 S (x) + k 04 , 
Ui(x) < kuSi(x) + ku, 
U r (x) < {k r i + k r2 S (x)}S r (x) + k r4l 



2 < r < r (1) , 

— — max' 



(2.7) 



and 



V (x) < k 03 Si(x) + k 05 , 
V r (x) < k r3 S p{r) (x) + k r5 , 

are satisfied, where 1 < p(r) < rmix for 1 < r < ril 



1 < r < r 



(2) 



(2.8) 



As a result of these assumptions, it is shown in [BL, Lemma 2.3 and Theo- 
rem 2.4] that, if r is such that 1 < r < rmL and if max{SV(xj^), S p ^(x^)} < 
C for some C, then there are constants C\ and C 2 , depending on C, r and T, 
such that 



sup S r (x t ) > C\ 

0<t<T 



< N- X C* 



(2.9) 



Semigroup assumptions 

In order to make sense of (II. 6ft . we need some assumptions about A. We 
assume that 

Aij > for all % ^ j > 0; < oo for all i > 0, (2.10) 
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and that, for some fi G R + + such that /i(m) > 1 for each m > 0, and for 
some w > 0, 

A T [i < wfx. (2.11) 
We then use [i to define the /x-norm 

:= $>(m)|ri on 7^ := {£ G ft: NCIU < °°L ( 2 - 12 ) 

m>0 

and, under these assumptions, the transition semigroup R is well defined 
[BL, Section 3], and 

^/i(i)i?y(t) < fi(j)e wt for all j and t. (2.13) 

Note that there may be many possible choices for /i, but that we also require 
that F is locally Lipschitz in the /z-norm, in order to ensure that (11. 6p has a 
/i-continuous solution: we assume that, for any z > 0, 

sup \\F{x)-F(y)\\J\\x-y\\ li < K(fi,F;z) < oo, (2.14) 
x +v- \\A\y-\\y\\v-< z 

and this should be borne in mind when choosing /i. We further assume that, 
for some (3 3 , /3 4 , 

fi{j) < "Of 3 and \A n \ < v{jf*. (2.15) 



Transition rate assumptions 

We need to ensure that the sum of the transition rates, even when weighted 
by largish powers of remains bounded. To ensure this, we assume that, 

(2) 

for some r large enough, there exist r% < Tmax, b > 1 and ki,k 2 > such 
that 

J>j(z)^|^|R7')r < {hS ri (x) + k 2 } b ; (2.16) 



JeJ 



this assumption is a specialized version of [BL, (2.25)]. In view of (12. 9p . this 
implies that, if max^^a;^), Sp{ri){ x i))} < C> then there are constants C\ 
and C 2 depending on C and T, such that 



sup 



0<t<T 



aj{x t i 



JeJ 



j>0 



< N- l C?. 



(2.17) 



We shall therefore assume that the initial condition needed for ( I2.17P is indeed 
satisfied: that, for some C < oo, 

max{S ri (<),S; (ri) (<)} < C. (2.18) 

It can be seen from the statement of Theorem 16.11 that the larger we can 
take r in (j2.16p . the sharper the approximation bound that we get in (II. 7p . 
in that ( can be taken smaller for a given value of the product (ro, resulting 
in larger values of b\ (£). 

Since it is immediate that 

it follows that, for any r, s > 0, 

J£Jm -HJm 

Je J i>o 

so that, if r + s < ro, (I2.17P implies that 

sup V {u{J)Ya.j(x? ) < dM- s , (2.19) 

except on an event of probability at most C 2 N~ 1 . 

Smoothness assumptions 

We need some smoothness conditions on the rates near the deterministic 
path x. First, for some 5 > 0, we assume that F has second order partial 
derivatives in the tube 

B(t,x,5) := {z E X: \\z -xJ M < 5 for some < s < t}, (2.20) 

where x solves (11.61) . and that, for any j, k, I, 

sup \D kl F j (z)\ < v jkl , (2.21) 

z£B(x,t,8) 

where the Vjki are such that 



$>(J>^ < K F2 n(k)n(l), (2.22) 
i>o 
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for some K F2 < oo. Note that (I2.22p is satisfied if 

Vjki < v J M^)M0> where < oo. (2.23) 

It is also true under the following condition: that, for each k, there exists 
N(k) C Z + with |iV(A;)| < n and ma.Xj £N ( k ) /i(j) < K Q fi{k) such that 

Vjki < Ki{fx(l)l {N ( k )y(j) + n(k)l{ N (j)}(j)}; v jk i = otherwise, 

(2.24) 

for suitable n ,K and Ki, all finite. The first derivative of F has already 
been assumed to be /x-Lipschitz in (12. 14ft ; with the assumption (12.211) . F 
becomes continuously /x-differentiable in the tube, so that, for some con- 
stant Kpi, 

sup V fx(j)\D k F j (x t )\ < K F1 fx(k) for all k. (2.25) 

We also assume that the individual transition rates aj are uniformly [i- 
Lipschitz in B(T,x,5), with 

sup \aj{z 1 )-a J {z 2 )\/\\z l -z 2 \\ II < K a {u{J)}^ (2.26) 

z u z 2 eB(T,x,5) 

for some K a , /3 5 > 0. This assumption, and those on the second derivatives 
of F, go beyond what is required for the law of large numbers in [BL]; the 
same is true of the assumptions (12.21) and (12.151) . 

Preliminary conclusions 

We now assume, in addition, that we can take 

r > 2(& + & + &) (2-27) 

in (12.161) . Then, under the assumptions of this section, it follows from [BL, 
Theorem 4.7], with ((j) := {v{j)} n \ that the following result holds: for each 
T > 0, there exist constants K^\k^ and such that, for all N large 
enough, if 

lh"-*„t<if«y^, (2.28) 

then 

P ^-^>K?^)<K ? ^. (2,9) 
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We shall from now on also assume that (I2.18P holds with x for x . Since 
then x can be represented as a limit of processes x M satisfying (j2.18p . because 
of (|2.29p . it follows in view of (I2.16P and (I2.17P that we also have 



sup V aj (x ( ) V|J J |Mj)} ro < C 1} (2.30) 
and therefore, as for (12.191) . for r + s < r , 



sup Y. {"(J)Yaj(x,) < C,M-. (2.31) 
a<t<T ^ 



When approximating by a deterministic path x, it is natural to choose 
their initial values to be close, as in (12. 281) . The impact of also assum- 
ing (I2.18P for the initial values of both paths is to specify how much closer 
the components need to be, whose indices j have large. 

Example 

In the model of Arrigoni (2003) presented in the introduction, we can 
take = j + 1, in which case /?i = 1 and /?2 = 2, the latter because of 
the migration transition. Calculation shows that (12.71) is satisfied for all r, 
as is (12. 8p also, with p(r) = 2r, so that we can take r^lx = Tmix = oo. 
Furthermore, (I2.16P is satisfied for any r , with r\ = r + 1. The quantities 
A and F are given by 

Ai = -{n + i(bi + d i + 'j)}; Aj^ = 2(^ + 7); Aj ;i+1 = ib h i > 1; 

with all other elements of A equal to zero, and, writing s(x) := X^>i i 2 ^ 
F\x) = P7(x i ~ 1 -x i )s(a;), i > 1; F°(x) = -p7x°s(x) + «, 

where we have used the fact that X^>o xJ = ^- Hence Assumption (I2.10p is 
immediate, and Assumption (12. lip holds for fj,(J) = j + 1 (so that {3 3 = 1), 
with w = maxj(6j — d{— 7 — /c) + (assuming the frj's and dj's to be such that this 
is finite). The value of 04 depends on the particular choice of the bi and di. 
For instance, the stochastic version of Ricker's (1954) model has both the bi 
and the di uniformly bounded, in which case we can take ^4 = 1. However, 
in the stochastic analogue of Verhulst's (1838) logistic model, the di grow 
linearly with i, and then one needs /3 4 = 2. 
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With the above choice of fi, F can easily be seen to be locally Lipschitz in 
the //-norm, with K(n, F; z) < Ap'yz. The partial derivatives of F are given 
by 

D k F\x) = pi k(x l - l -x l )+ P1 s(x){l {k} (i-l)-l [k} (i)}- 
D kl F\x) = p7{A;[l w (i-l)-l w (z)] + i[l {fe} (z-l)-l {fc} (z)]}, 

for any i, k, I > (we take x~ l = 0). From this, it follows (using the elemen- 
tary bound j + 1 > 2j in j > 0) that we can take K FX = 6p7sup 0<t<T 
in (12.251) . and that (12.241) is satisfied with n = 2, K = 2 and Ki = p7, so 
that (12.221) is also satisfied (one can in fact take Kp 2 = 4p7). Finally, ( I2.26P 
is satisfied, with (3^ = if the di are uniformly bounded, and with f3§ = 1 if 
they grow linearly, and with K a of the form K'{1 + sup 0<t<T \\x t 

3 Controlling r] N 

From now on, we assume that all the assumptions of Section [2] are in force. 
We first show that the effect of the perturbation r] N is negligible. For rj^ , 
from ( I1.13p . we need to consider the difference 

/ R(t - s){F(xf) - F(x.) - DF(x s )[xf - x s }} ds. 
Jo 

We note first that, if \\h\\,, < 5 for 5 as in Condition fl2~^Uj) . then, from (12"^2"|) . 
\\F(xf)-F(x s )-DF(x s )[h}\\, < ^MJ)EEI^I^ ^ K F2 \\h\\, 2 . 

j>0 k>0 l>0 

Hence, from f!2.29|) and from (I2.13p . for all N large enough to ensure that 
KlpN-V 2 y/\ogN < 5, we have 

sup ||??f|^ < Kp 2 t{K < >p} 2 e wt N~ 1 l' 2 log N, (3.1) 

0<s<t 

for all < t < T, except on a set of probability at most K^N" 1 log N. 

4 Discrete to diffusion 

We now show that N l l 2 m N is close in the p-norm to the diffusion W, given 
by 

W t := ^JWj(Aj(t)), 
13 



as in (j!.14p . We first need to show that this W indeed has paths in IZ^. For 
this, it is enough to show that 

5>(i)El ji IW4^))l < oo (4.1) 

i>o Jej 

for all t. 

We begin by noting that, using the reflection principle, if B is standard 
Brownian motion, then there exists a constant 7 < 00 such that, for all a > 0, 

P[sup \B(x)\>arf] < e 1 "^ . (4.2) 

0<x<l 

Thus, from (I4.2p . for any C > 1 and p, T > 0, we have 



\Wj(Aj(t))\ = d Aj(Ty/ 2 \B(Aj(t)/Aj(T))\ < Aj(T)V 2 Wplog(Cv(J)), 

for all < t < T, except on a set of probability at most e{CV(J)}~ p . Hence 
it follows that 



\Wj(Aj(t))\ < AjiTY^y/plogiCviJ)) (4.3) 
for all < t < T and for all J G J , except on a set of probability at most 

For x,y > 0, one has y/x + y < (1 + + ^/y). Substituting from 

( 14. 3 p into (14. ip shows that || W*||^ < 00 a.s. for all < t < T, provided that 

j>o Jej 

since C is arbitrary. However, by (I2.15P and recalling the definition of J*, we 
have, for any e > 0, 

Emj)Ei j >( j ) £ ^) 1/2 

j>o Jej 

< \J 3 njy M2 - £ ^(j) ft+2£+/32/2 A J (t) 1 / 2 



Jej j>0 



V2 / T \ 1/2^ 



< J J PT K^ 2 ^ I E / z/(J)^ +2 ^+^ aj ( Xs ) rfs j , (4.4) 



^JeJ / \JeJ 
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and both sums in the final expression are finite, by (12.41) and (12.301) . provided 
that 02 + 2/?3 < r and that e is small enough. 

Having established that W indeed has paths in TZ^, we now need to show 
that it is close to N l / 2 m N in the /x-norm, if the Brownian motions Wj are 
suitably chosen. The relationship between N l l 2 m N and W arises because 
can be represented in the form 

iVmf := n\x?-x$- 



= ^{P^NA^-NA"®}, (4.5) 
Jej 

where Aj(t) := J * aj(x^) ds, and the P/s are independent Poisson pro- 
cesses. Now {N~ 1 / 2 (Pj(Nt) — Nt), t > 0} can be well approximated by 
a Brownian motion, and Aj(t) is close to Aj(t) := f aj(x s ) ds, by f !2.26|) 
and (l^29|) . 

We thus wish to show that the Wj can be chosen in such a way that 



sup 

o<t<r *rrz 
— j>o 



J2j J Zy(A»(t))-J2j 3 Wj(Aj(t)) 

Jej Jej 



(4.6) 



is small, where we define 

Zf(s) := N' 1 / 2 {P J (Ns)-Ns}. (4.7) 

For use in the next section, we prove somewhat more: that, under ap- 
propriate conditions, we can replace fi(j) in ( I4.6P by the larger quantity 
'■= { Z7 (j)} /33+/34 ; an d still obtain something that is small. To do so, we 
begin by bounding the sum by T x (t) + T 2 (t) + T 3 (t), where 

T i(t) ■= E Y.\ Jj \^\ z ^ AN j^))- w Mj(t))\\ 

JeJ M j>o 

m) := E El JJ >*0')l^(^(*))i; (4-8) 
■HJ M i>o 

T *(t) ■= E Ei JJ >*0')i^(^(*))i- 

Here, M is to be chosen later as N^, for some suitable small £ > 0. 
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We begin with T 2 (t), which we deal with by showing that, for suitable 
choice of M, NJ2mj m Aj (T) is small. Indeed, 

N E = N f E "'CO**, 

which is bounded by using (I2.19P with r = and s < r , together with (12. 17ft ; 
the quantity is of order NM~ S for any s < Tq, except on an event of proba- 
bility of order 0(iV _1 ). Thus, except on an event with probability of order 
0(N~ 1 + NM~ r °), N~ l l 2 Pj(NA^{T)) = for all J (£ J M . Furthermore, the 
contribution from the compensators is bounded by 

N 1/2 E <j)\J j \A N j{T) < N^T sup E ^* +/ V"M*T) 

0(N 1/2 M~ 3 '), (4.9) 



JiJ M °- t - T HJm 



by (I2.19p . if /3 3 + /3 4 + s' < r , except on an event with probability of or- 
der 0(iV _1 ). Recalling (12.271) . this proves that, for any s' < r — (3 3 — /3 4 , 

sup T 2 (t) = 0{N 1/2 M~ s '), (4.10) 

0<t<T 

except on an event of probability of order 0(iV -1 + NM~ r °). 
For T 3 (t), we use (14. 2 p to give 



sup \Wj(Aj(t))\> {Aj^Y^ploguiJ) < e{u(J)}-P, (4.11) 

0<t<T J 



P 

for any p > 0. Hence, for any p > (3 2 , it follows that |W/j(v4j(t))| < 
{A j(T)} 1 / 2 7A/plog u(J) for all J ^ Jm and for all < t < T, except on an 
event £"3 of probability of order 0(M~ p+l32 ), from (12. 4p . But then, except 
on E 3 , for all < t < T, 

E Ei J >*^i^(^))i ^ EEi j >«w^( t )}^vpW) 

< K £ ^> E JMJf' M£ {MT)} l/ \ 
J$Jm 

for any £ > 0, with suitable choice of K £ . But now, for any r > 0, 
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n 1/2 , >> 1/2 

HJm J \HJm 

< {^-^-^^{^^-^^l^^M-^} 172 , (4.12) 

by (J23D and f l2~3T|) . so long as r > f3 3 + /3 4 + e + /V 2 and 2(r + r') < r . 
Hence, if r > /3 2 + 2(J3 3 + /3 4 ), then for any r' < (r - /3 2 )/2 - (/3 3 + /3 4 ) we 
have 

p- 1/2 sup T 3 (t) = 0(M~ r '), (4.13) 

0<t<T 

with an implied constant uniform for all p > 02, except on an event with 
probability of order 0(M~ P+I32 ). 

So far, the bounds have been achieved without any specific choice of the 
Brownian motions Wj, but, for Ti(t), we need to be more precise. We treat 
each J separately, since the underlying Poisson processes Pj are independent, 
and match the centred and normalized Poisson process to a Brownian 
motion Wj using the KMT construction. We need only to do this over a 
limited time interval, since, from (I2.17p . 



sup E a ^( x f) ^ Cij 



except on an event E of probability of order 0(iV _1 ), so that, off E , 

A j (T) < TC X for all J. (4.14) 

We use Komlos, Major & Tusnady (1975, Theorem 1 (ii)), together with (14. 2 p 
to interpolate between integer time points, applied to the centred unit rate 
Poisson process. This implies that, for any p > 0, we can choose Wj in such 
a way that 

sup {N~ 1/2 (Pj(Nt) -Nt) - Wj(t)} < k p N- 1/2 \ogN, 

0<t<TCi 

for a constant k p , except on an event Ej p of probability of order 0(N~ P ). 
Thus the same bound holds for all J G Jm except on an event E p of proba- 
bility of order 0(M^ 2 N~ p ). Hence, except on E UE p , an event of probability 
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of order C^iV" 1 + M^N~ P ), we have 

T n(t) ■= £ £ i^'l \Z?(Aj®) ~ Wj{A»(t))\ 
JeJ M j>o 

JeJ M 

= 0(M Ms+ ^N- 1/2 \ogN), (4.15) 
for all < t < T. It thus remains to bound 

Tu(t) := E Y,\ jj \^)\ W M N j{t))-Wj{Aj{t))\- (4-16) 

First, note that, by f[2T26]) . 

sup < T^ Q {z/(J)} & sup Hxf-xtlU (4.17) 

0<t<T 0<t<T 

and that 

sup \\xF-xt\U < K ( T } N- 1/2 ^/\ogN, (4.18) 

0<t<T 

by (I2.29p . except on an event with probability of order 0(N~ l log N). Fur- 
thermore, by (I2.30p . Aj(T) < TC\ for all J. Then, for a Brownian motion W 
and for < 8 < 1, by a standard argument based on (I4.2p . 



< {A5~ 1 + 2)e5 r 



sup |W(u + s) — W(«)| > 375 1/2 v / rlog(l7^) 

0<u<A,\s\<8 

(4.19) 

for any r > 0, with 7 chosen as for (14.21) . Hence, taking W = Wj and, in 
view of f HTTj) and (gjgj> , taking 



5 = 8j = TK a K^{u(J)}^N-^ 2 ^hg~N 

and A = TC\ in (I4.19p . it follows that, for any s, r' > 0, there is a K r i £ < 00 
such that 

sup \Wj(A N j(t))-Wj(t)\ < K r , £ M^ /2 N- 1/i+£ for all J e J M , 

0<t<T 
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except on an event of probability of order O^N^ 1 log N+M h+2 ^{M 2 ^ /N} r ') 
Off the exceptional event, we have 



sup T 12 (t) < J* V u^K^M^N' 1 ^ 

= 0(M /32+/33+/34+/35 / 2 AT- 1 / 4+e ), (4.20) 



for any e > 0, and the exceptional event can be made to have probability of 
order 0(N~ l logiV) by choosing r' large enough, provided that M is bounded 
by a small enough power of N. 

Combining ( jOUj) . (ETT3"]) . (ETT5]) and (ET20]) . and choosing M = iV^, we 
see that we have no useful bound unless (r > 1 (because of the exceptional 
event in (00) ) and C < l/{4(/3 2 + /3 3 + /3 4 ) + 2/3 5 } (in view of (OP]) ), so 
that r > 4(/3 2 + /?3 + ^4) + 2/3 5 is a minimal requirement. Note that this 
assumption on r is more restrictive than that assumed in Section [2J The 
error bound in (I4.15P is always smaller than that in (14.20) . and with (ro > 1, 
the error bound in f |4.10j) is smaller than that in (I4.13p . This translates 
into the following conclusion: if r > 4(/3 2 + 03 + 0±) + 2/3 5 , then for any 
l/r < C < l/{4(/? 2 + 03 + 0a) + 20 5 } we have 

sup HiV^mf - V^IU 

0<t<T 



< sup VMj)} 

0<i<T — ' 



j>0 



Y J J ] Z^A N J (t))-Y J J 3 W J {A J (t)) 

= 0(N~ bl ) (4.21) 
except on an event of probability of order 0(N~ b ' 2 ), for any 

h < 61(C) 

:= min{i-C(/3 2 + /33 + /34 + /3 5 /2),K(^-/32-2(/33 + /34))}; (4.22) 
b 2 < b 2 (() ■= min{Cr - 1, 1}. 



5 The existence of W, and its approximation 

The next step in the argument is to show that the process W in (II. 15) . 
related to W exactly as m N is related to m N through (ll.lOp . is well defined, 
and that it is indeed the limiting analogue of the process N 1 ' 2 ?^ . For its 
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existence, recalling (11 . 15)) . it is enough to show that R(t — s)AW s exists for 
each s,t, and belongs to TZ^. For this, it is enough to show that 

5>(o Rad - *) E 1^*1 E Jk\wj(Aj( S ))\ 

i>0 j>0 k>0 Jej 

is a.s. bounded. Now, in view of (12. lip . (12. 13)) and (12. 15)) and recalling that 
the off-diagonal entries of A are non-negative, this will be the case if we can 
bound 

J2»(k)(i + \A kk \)J2\ J k\\Wj(Aj( s ))\ < j*Y,HJ)Y* +Pa \wMj{s))\ 

fc>o Jej Jej 

(5.1) 

uniformly in s. Now, once again from (I4.2j) . for any C, r > 0, 



sup 

0<t<T 



\Wj(Aj(t))\ > {A J (T)} 1 / 2 7v /rlog(Cz/(J)) 



< e{Cu(J)}- r , 



so that, in view of (12.41) . taking any r > 02, there is a (random) C such that 



sup 

0<<<T 



\Wj(Aj(t))\ < 7 VF{A J (T)} 1 /Vlog(Ci/(J)) 



a.s. for all J. But now, returning to ( 15. ip . we just need to show that the 
quantity 

E W J ) V 3+ ^ +£ {Aj{T) } 1 ' 2 

Jej 

is finite for some e > 0, and this is achieved as in (I4.12p . if r > /? 2 +2(/3 3 +/3 4 ). 

To show that W is a good approximation to we begin with the 

result proved in (14.211) above, that, except on an event of probability of or- 
der 0(N- b2 ), sup < t < r WN^mf - W t \\n = 0(N~ bl ) for any b x < h(C),b 2 < 
62(C)- This quantity is one element of WN^-^m^ — W t \\^; the other is 

/ R(t - s)A(N 1 / 2 m^ - W.) ds. 
Arguing much as in the previous paragraph, we need to bound 
sup ^2MJ)V 3+ ^\Wj(Aj(t)) - Z»{A N j{t))\- 



o<t<T JeJ 
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But this is exactly what we achieved in ( 14.2 ip . Hence, for ( such that l/r < 
C < l/{4(/3 2 + & + &) + 2/3 5 } and for any b x < MO, h < 63(C), 

sup WWt-N^m^W, = 0{N~ h ), (5.2) 

0<t<T 

except on an event of probability of order 0(N~ b2 ). 



6 The final approximation 

The final step in the argument is to compare the solution U to (I1.12p with 
the solution Y to (11.161) . Both satisfy the general equation 

Z t = R(t)Z + [ R(t-s)DF(x s )[Z s ] + z t , (6.1) 
Jo 

but with different initial conditions Z and forcing functions z; and their 
difference Y — U N also satisfies (I6.ip . with initial conditions and forcing 
functions subtracted. Now, for U N , the forcing function rj N + N 1 ^ 2 m N is 
close to the forcing function W for Y, because of (13.11) and (I5.2p . and we 
shall assume that Uq and Yq are also close to one another, so that both 
differences are small. We now show that this implies that the difference 
between Y and U N is also small. 

First, the assumption (I2.25P implies that, for w 6 TZ^ \\DF{x s )[w]\\^ < 
^fiIIHIai- ^ i s then immediate from (I2.13P that 



R(t - s)DF(x s )[Z s ] ds 



< K m / e w{ ~ t - 8 HzA\„ds 



Fl I e W^sWn' 



ft 

< K F1 e wt I WZ.W^ds, 





and that ||i2(t)Z ||^ < ||Z || At e w ' t . Hence, for < t < T, it follows that 



H^tlU < \ \\Z \\^e wl + sup \\z s \\Ae ct } (6.2) 

L 0<s<T } 

with C = K F ie wT . Applying ( 16. 2 p to Z = Y — U N , and using the bounds 
in ( 13. ip and ( 15. 2p . it follows that, except on an event of probability of order 
0(iV- b2 ), 

sup \\Y t -U t N \\^ = 0(N- bl ), 

0<t<T 
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where U N := N l ' 2 (x N - x) and Y is the solution to (jl.lGp . provided that 
|| Yo — Uq \\p = 0(N~ bl ) also. This proves the main theorem of the paper: 

Theorem 6.1 Suppose that the assumptions of Section^ are satisfied, and 
that we can take r > 4(/3 2 + P3 + 04) + 2/3 5 in (12.161) . For any ( such that 
l/r < C < l/{4(/3 2 + Ps + Pa) + 2p 5 }, define 

6i(C) := min{i-C(/3 2 + /33 + /34 + /35/2),|C(ro-/32-2(/3 3 + /34))}, 

andb 2 (C) ■= min{Cr - 1, 1}. Suppose that \\Y - U^W^ < KN~ bl ^. Then, 
for any bi < &i(C); ^2 < ^(C); we can construct copies ofY and U N on the 
same probability space, in such a way that 

sup ||y t -0f||„ = 0(N- b "), 

0<t<T 

except on an event whose probability is of order 0(N~ b2 ). 

So, for example, in the model of Arrigoni (2003), we can take r as big 
as we wish, and then (ro = 2, allowing b\ = 1/4 — e and b 2 = 1 — e for 
any e > 0. However, these rates can only be attained for correspondingly 
well controlled initial conditions: in addition to ( 12. 28ft , it is necessary to 
ensure that (12.181) is satisfied, so that S2( ro +i)( x o) < C for some C < oo and 
for all N, and that 5 , 2 ( ro+ i)(x ) < C also. For stochastic logistic dynamics 
within the patches, with bi = b and di = d + ci, we need to take ro to exceed 
4(/?2 + 03 + Pi) + 2/3 5 = 22 to yield an error bound that converges to zero 
with N, and thus require the initial conditions to have uniformly bounded 
(46 + S)-th moments for some 5 > 0. 
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